Heating mechanism affects equipartition in a binary granular system 
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Two species of particles in a binary granular system typically do not have the same mean kinetic 
energy, in contrast to the equipartition of energy required in equilibrium. We investigate the role 
of the heating mechanism in determining the extent of non-equipartition of kinetic energy. In most 
experiments, different species are unequally heated at the boundaries. We show by event-driven 
simulations that differential boundary heating affects non-equipartition even in the bulk of the 
system. This conclusion is fortified by studying numerical and solvable stochastic models without 
spatial degrees of freedom. In both cases, even in the limit where heating events are rare compared 
to collisions, the effect of the heating mechanism persists. 



An issue of broad generality for noncquilibrium physics 
is the assignment of intensive temperature variables to 
steady states of driven systems. A widely-studied exam- 
ple of practical importance is a dilute system of macro- 
scopic grains. A granular gas of this kind is different 
from an equilibrium gas in that a continuous supply of 
energy from an external source is necessary to maintain 
a steady state and balance dissipation due to inelastic 
collisions between particles. The energy source deter- 
mines typical scales of the dynamics in the system, but is 
presumed not to affect the microscopic constitutive laws 
such as the equation of state or expressions for trans- 
port coefficients In this Letter we study an externally 
nuidised granular gas and find that the details of the 
energising mechanism affect relations between intensive 
quantities even in the bulk of the system. 

In a typical real-world situation, energy is delivered at 
the boundaries of a granular system by vibration, shear or 
other mechanical means. In equilibrium, a gas placed in 
contact with a heat bath acquires a uniform temperature 
and density. However, in the case of a granular system, 
there are gradients in density and particle motions as a 
function of distance from the energizing boundary. For 
dilute gases of inelastic particles there has been consider- 
able progress 0]in describing this inhomogeneous steady 
state in terms of density and "temperature" fields, where 
the temperature is a purely kinetic construct, defined as 
the mean kinetic energy per degree of freedom, by anal- 
ogy to kinetic theory of a molecular gas. 

An important conceptual tool in validating the use 
of nonequilibrium temperatures involves thermal contact 
between systems @i H|. One can ask whether a candi- 
date definition for the temperature governs the direction 
of energy flow and whether the temperatures equalise or 
evolve to bear a fixed relationship. In this spirit, we 
consider a granular system with two species of particles. 
There is no requirement that the species acquire the same 
kinetic temperature; indeed, experiments @, @, 0] docu- 
ment a violation of equipartition. It is observed that 
two species a and b acquire temperatures whose ratio 
7 = T a /Tb is affected strongly by the ratio of particle 
masses, m a /nib, and weakly by their inelasticity. The 



temperature ratio 7 is observed [6| to be only weakly 
affected by the composition - the number density and 
stoichiometry - of the mixture. These observations have 
been reproduced and more comprehensively explored by 
event-driven simulations [1, 0, [1(3] of vibration-fluidised 
systems. 

This raises the question of whether the extent of depar- 
ture from equipartition can be predicted, that is, can the 
temperature ratio 7, be expressed in terms of specified 
properties of the two types of particles in the mixture. 
If so, it would not be necessary to introduce indepen- 
dent temperature fields to describe each species in the 
mixture. Following early work ll| that addresses this 
question within the framework of kinetic theory, there 
have been several theoretical treatments of this issue in 



spatially homogeneous models. Garzo and Dufty [12 1 
studied theoretically the homogeneous cooling of a bi- 
nary granular mixture and found that cooling rates of the 
two species are the same while their temperatures remain 
different, thus maintaining a time-independent 7. These 
predictions have been compared against MD simulations 
[l3(. Likewise, kinetic theory treatments of homoge- 
neously heated steady-states of mixtures have been suc- 
cessfully compared to event-driven simulations in which 
particles get occasional velocity boosts from an external 
noise source 14 , [H| . Other treatments of the homoge- 



neously heated case have been advanced by Langevin- 
like models [lfj} and by mean-field solutions 17 1 of the 
Boltzmann equation. All of these theoretical results have 
been compared to the 7 found in spatially inhomogeneous 
vibration-fluidised steady states [1, @, 0, S, EJ and 
while some qualitative agreement has been found, it is 
appropriate to be cautious in comparing these situations. 

A further, and often-neglected ingredient in this prob- 
lem is that in a typical experiment, the driving mecha- 
nism couples differently to different species of particles. 
For instance, in a vibration-fluidized system [H, la, 0] , the 
vibrating boundaries may be treated as infinitely mas- 
sive, and therefore, impart a characteristic velocity scale 
(rather than an energy scale) to the particles. If the par- 
ticles in the system are of different mass, then they are 
differentially heated at the boundary. This is not an is- 
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sue in an equilibrium system, where equipartition holds 
even if the two species are coupled differently to a heat 
bath. In the granular systems under consideration, the 
temperature ratio 7 is experimentally observed [y] to be 
close to the mass ratio of the two species of particles 
near the boundary, and then attain an apparently con- 
stant value after a few mean free paths. The question 
remains whether this apparently constant bulk value de- 
pends on the boundary condition. A similar issue has 
previously been commented on in the context of a ho- 
mogeneously heated system with both species coupled to 
different heat baths [15] and in comparing homogeneously 
heated systems to boundary-driven systems 0]. 

In this Letter, we address by molecular dynamics simu- 
lations the effect of this differential boundary heating on 
the extent of nonequipartition of energy in the bulk of 
the system for a binary granular mixture. We find that 
the level of differential heating affects nonequipartition 
in the bulk, and that the effect of the boundary is never 
forgotten. These results are strengthened by considera- 
tion of a numerical simulation model [l8| without spatial 
degrees of freedom, and of a stochastic modelflij. 

In our event-driven simulations, two types of smooth 
discs (labelled a and 6), both of diameter d, but with 
different masses, are fluidized inside a rectangle of di- 
mensions 48 x 32d. These simulation parameters are 
chosen to mimic an experimental geometry previously 
studied by us @. For the simulations reported here, the 
ratio of particle masses rrib/m a is 5, the restitution co- 
efficient of both species is 0.93 and an equal number of 
both species are used with the total number of particles 
2N=200 corresponding to area fraction 10.2%. We set 
.9 = and employ periodic boundary conditions in the 
horizontal (x) direction. In order to isolate the effects 
we are concerned with, we choose to heat the particles 
using fixed "thermal" walls rather than with the addi- 
tional space- and time-dependence introduced by oscil- 
lating walls. Following a collision with the top or bottom 
boundary, a particle is reflected with normal and tan- 
gential velocity components independently drawn from 
gaussian distributions. The width of the gaussian is cho- 
sen to be < v% >=< Vy >= 2E a ^/m a ^ for the two types 
of particles, thus allowing us to control the differential 
heating at the boundary. 

In FigQ]we show three choices of differential heating: 
Eb/E a = I where equal amounts of energy are fed to 
two types of particles at the boundaries, Eb/E a = 5 
where the energy fed is proportional to the mass ra- 
tio as in the vibrational experiments, and Eb/E a = 25 
which even more disproportionately emphasizes the heav- 
ier particles. FigQ] displays (a) the average tempera- 
ture ratio 7 = T a (z)/Tb(z) and (b) number density ratio 
n a (z) I 'rib(z) profiles along the direction perpendicular to 
the top and bottom boundaries. Comparing any single 
pair of curves - say the case Eb/E a = 5 - with experi- 
mental data Q we find similar profiles of 7 and number 
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FIG. 1: (a) Temperature ratio profile and (b) number density 
ratio profile for two species of particles along the vertical (z) 
direction. The particles differ only in their masses with nib = 
5m a . There are 100 particles of each species in a rectangle 
of size 48 x 32d, corresponding to an area fraction of 10.2%. 
Three levels of differential heating are shown: E a = Eb (dotted 
fine), E a = i?b/5(solid line), E a = -E(,/25 (dashed line). 



density ratio, although the energy-supplying boundaries 
are fixed and gaussian. In each of the three cases, the 
value of 7 at the boundaries is determined by the heating 
mechanism, as expected, but changes to a nearly constant 
value in the bulk. However, taken together, the three 
cases show clearly that the bulk value of 7 depends on 
the choice of differential heating at the boundary. Thus 
shows that the value of 7 is not an intrinsic property 
of the pair of particle species and that the effect of the 
boundary permeates the bulk of the system. 

We further check whether the persistence of the bound- 
ary condition into the bulk is a finite size effect by dou- 
bling the height of the system to 64c? while holding fixed 
the total area fraction. The three heating schemes shown 
in Fig[T] are once again employed, and the temperature 
ratio and number density ratio as a function of z are 
shown in Fig[2] Clearly, the temperature ratio 7 is still 
affected by the value of Eb /E a specified at the boundary. 

It is straightforward to simulate substantially taller 
systems, however, the increasing segregation of the heav- 
ier particles to the middle of the system presents a new 
complication. We thus address the issue of whether the 
data in Fig[T] and [2] are still affected by finite-size limi- 
tations by considering the number of collisions q suffered 
by a particle since it has collided with a boundary. Each 
time a particle i collides with another particle, qi is in- 
cremented by 1, and after it collides with either heating 
boundary, qi is set to zero. The probability distribution 
of the vertical distance I of a particle from the wall it 
last collided with, is shown in the inset of Fig[3]for var- 
ious values of q. As expected, P(l) broadens and shifts 
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FIG. 2: Vertical temperature ratio profile (a) and density ra- 
tio profile (b) in a system of height 64d, double that of the sys- 
tem in FigJT] with the area fraction maintained at 10.2% The 
lower panel shows greater segregation of the massive species 
to the centre of the system than in Fig[T] 
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FIG. 3: The temperature ratio, 7, as a function of collision 
number, q, under three levels of differential heating. The 
corresponding dotted curves sample only those trajectories 
that start at one heating wall and terminate at the other. 
The inset plots the distribution of the distance from the last 
heating boundary I for different ranges of q. 



monotonically to larger I, for larger q, saturating at the 
half-height of the system 16d for large q. Next, statistics 
for the kinetic energy of the two types of particles as a 
function of q are sampled at fixed time intervals during a 
long simulation time and the temperature ratio 7 is ob- 
tained as a function of q. As shown in Figl3l each curve 
of 7(g) reaches a plateau after a small number of colli- 
sions. This is consistent with MD simulations of bidis- 
perse systems [Hi] where it was found that after about 
10 collisions the temperature ratio 7 reached a steady- 
state value from arbitrary initial conditions. However, as 
can be seen in FigO the asymptotic value of 7 at large 
q remains a function of Ei>/E a . Therefore, details of 
the driving mechanism at the boundary appear not to be 
erased by interparticle collisions. To ensure that this con- 
clusion is not dominated by particles that suffered many 
collisions but happened to remain close to a boundary, 
we also plot in FigJ3]7 for the subset of particle trajec- 
tories that start at one heating wall and terminate at 
the other. The plateau values of 7 of these trajectories 
coincide with those from the entire set of particles. 

It may still be argued that a finite-size effect persists 
since particles with large q suffer collisions with particles 
that have recently collided with the boundary. In order 
to further confirm the persistent effect of the differential 
heating mechanism we turn to a numerical model with- 
out spatial degrees of freedom, based on that introduced 
by van Zon and MacKintosh||lj| to simulate a monodis- 
perse system. A system of 2N(N=100) inelastic particles 
is initialized with gaussian-distributcd initial velocities, 
and in each time step, C pairs of particles are randomly 
selected to collide in a 2-dimensional inelastic collision, 
and H particles are selected for heating. The impact 



parameters for the collisions are chosen from a uniform 
random distribution. For each heating event, velocities 
are selected randomly from Gaussian distributions with 
average energy E a and depending on the species of 
particle selected. Once again, the three cases of differen- 
tial heating used in the spatial event-driven simulations 
can also be implemented here. For each of these three 
cases, statistics for the temperature ratio 7 are accrued 
for many values of q avg = 2C/H, the ratio of the fre- 
quency of collisions to heating events. (The large q avg 
limit corresponds to a boundary driven system in which 
the inter-particle collisions are much more frequent than 
heating events, and corresponds to large q in FigJ31 with- 
out the complications introduced by spatial gradients). 
As shown by the three sets of data symbols in FigJU the 
temperature ratio 7 as a function of q avg is different for 
the three heating mechanisms even for large q aV g, where 
the heating events are rare relative to the collisions that 
one might intuitively expect to erase details of the heat- 
ing mechanism. 



The numerical model discussed above can be described 
in terms of an exactly solvable stochastic mode l 11911 . We 
have adapted the model of Shokef and Levine[19( to in- 
corporate equal numbers of two types of particles a and 
b. The time evolution of the energy of a particle 

i of either species during an infinitesimal time step dt is 
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expressed by a stochastic equation, 
value: 



1.0E 



E l {t + dt) 



probability: 
1 - Pdt, 
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^ or El, 



(1) 



(1 - /)Pdt, 
fPdt, 



where E 1 ^ (t) is the instantaneous energy of the particle 
at time t, P is the interaction rate per particle per 
unit time, / is the fraction selected for external heating, 
Xij—mi/rrij, r) is the restitution coefficient, E l ab is the 
external energy given to particle i for a heating event 
with an input energy E % a or E\ depending on the species 
of particle i. Therefore, the first line on the RHS of Eq[T] 
corresponds to particle i not undergoing interaction, the 
second line corresponds to particle i being selected for 
collision with another particle j (averaged over collisions 
with species a and b) and the third line corresponds to 
particle i being selected for heating. For the second line, 
we use a simple one-dimensional collision model with un- 
corrected pre-collisional velocities. In steady state, the 
ensemble averages obey < E l (t + dt) >=< E l {t) >. By 
solving this stochastic equation, we determine the de- 
pendence of the temperature ratio 7 on q avg as shown in 
FiglU To make comparison to the numerical simulations 
of the random heating model, we make the correspon- 
dence / = H/(H + 2(7) = 1/(1 + q aV g)- Once again, 
three different values of Eb/E a are shown. Despite the 
simplified treatment of collisions within our adaptation 
of the Shokef-Levine model, we obtain reasonable agree- 
ment between the numerical model and this stochastic 
model. Qualitatively, the dependence on Eb/E a follows 
the same trend, and the smooth curves in FigJJ] for the 
three heating conditions do not converge at large q aV g- 

In summary, we have demonstrated by three indepen- 
dent means, the persistent effect of the boundary heat- 
ing mechanism on the extent of nonequipartition in a 
binary granular system. The specific context for study- 
ing the effect of differential heating come from studies of 
vibration-fluidized granular systems, however, this type 
of differential heating should be generic to other forms 
of driving: in shear, for example, particles will be dif- 
ferentially excited depending on their frictional proper- 
ties, shape, or size relative to the features of the shear- 
ing surface. Similar effects may be anticipated even in 
monodisperse systems where the equipartition between 
rotational and translational degrees of freedom may not 
be determined purely by particle properties but by the 
degree to which each degree of freedom is pumped by the 
driving mechanism. More generally, details of the bound- 
ary heating mechanism can not be ignored in describing 
inelastic gases, leading to concerns about quantitative 
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FIG. 4: The temperature ratio 7 as a function of the aver- 
age particle-particle collision times q avg - Numerical model 
results are expressed in lines with symbols, E a = £4, (open 
circle), E a — i5 b /5(cross), E a — i? b /25(solid square). The 
smooth curves correspond to results derived from the stochas- 
tic model, E a = (dotted line), E a — E /5 (solid line), 
E a = -Eft/25 (dashed line). The corresponding asymptotic val- 



ues as q avg 



00 for these curves are 0.926, 0.798 and 0.755. 



comparisons between theories of homogeneously excited 
granular systems and boundary-driven experiments. 
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